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It is demonstrated that an extra nodal point and a domain with anomalous dispersion are present 
in the spectrum of charge carriers in graphene, when the concentration of weakly bound impurities 
exceeds a certain critical value, determined by the spatial overlap of individual impurity states. 
Corresponding spectrum rearrangement is shown to be of the cross type. 
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I. INTRODUCTION 

Chemical functionalization of graphene is currently a 
rapidly developing field of researchJ^ By depositing var- 
ious atoms, molecules and chemically active groups on 
graphene, it is possible to tune up its properties ac- 
cording to specific technical applications. Undoubtedly, 
transport properties of graphene arc decisive for its fu- 
ture in microelectronic devices of the new generation. 
With respect to the transport properties, adatoms should 
be viewed, first of all, as a sort of defects, which can 
significantly alter the spectrum of charge carriers and 
lead to the localization of states. It was demonstrated 
experimentally that ion bombardment 3 -, deposition of 
atomic hydrogen^ or fluorine^ opens a mobility gap in 
the electron spectrum of graphene and triggers the metal- 
insulator transition, just as it happens in conventional 
systems with massive carriers. The opening of the mobil- 
ity gap in graphene is attributed to the presence of impu- 
rity resonance states in the Dirac point vicinity. More- 
over, it became evident that some amount of resonant 
defects is always present in graphene obtained by the 
micromechanical cleavage process^ It follows from theo- 
retical studies of generic models for adsorbed atoms that 
such defects should considerably modify the electron den- 
sity of states and affect the conductivity of graphenej 7 - - — 
However, only those cases were examined, in which im- 
purity levels of adatoms were strongly bound to the host. 

Lifshitz defects (vacancies including), which do noth- 
ing to the system but change effective potentials on lat- 
tice sites occupied by them, bear a lot of similarity to 
strongly bound impurity centers. 12 Albeit such defects 
are capable in producing resonance states in the Dirac 
point vicinity, the corresponding resonances are of a re- 
duced integral intensity, and the subsequent spectrum 
rearrangement has a diffuse, anomalous character, which 
is inherent in low-dimensional systemsi 13 ' 14 This type of 
the spectrum rearrangement is characterized by the pres- 
ence of two non-overlapping dispersion branches in the 
spectrum, which are separated by a mobility gap. At 



that, the density of states inside this gap can consider- 
ably exceed the one in the hostJ^ 

However, there is a different opportunity: spectra of 
elementary excitations in tree-dimensional systems (and 
in low-dimensional ones as well) can undergo a rear- 
rangement of the coherent (cross) type with increasing 
the amount of impurities^ The renormalized disper- 
sion resulting from the spectrum rearrangement of this 
type looks similar to the hybridization between the host 
branch and the dispersionless branch that corresponds 
to the impurity state energy. Namely, two overlapping 
branches develop in the spectrum of a disordered system. 
As a consequence, two different energies of quasiparti- 
cles correspond to each wave vector. These dispersion 
branches are separated by a true gap, which is gradu- 
ally broadening with increasing the impurity concentra- 
tion. The cross-type spectrum rearrangement is usually 
accompanied by strong single-impurity resonances, which 
appear in low-dimensional systems only for weakly bound 
impurities. Thus, we examine below possible changes in 
the electron spectrum of graphene, which can occur under 
an increase in the impurity concentration, when impurity 
levels are weakly bound to the host. 



II. IMPURITY MODEL 

When dealing with 7r-electron bands of graphene 
within the tight-binding approximation, it is usually suf- 
ficient to take into consideration only matrix elements 
between nearest neighbors in the honeycomb latticed 



H — Hq + Hi, 
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where H is the dimensionless host Hamiltonian, H im 
is the impurity perturbation, vector n runs over lattice 
cells, indices a and j3 enumerate sublattices, c^ Q and c na 
are the creation and annihilation Fermi operators at the 
corresponding lattice site, and the prime at the sum sign 
indicates that na ^ mj3. For the sake of simplicity we 
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choose the energy unit in such a way that the transfer 



integral t = 1/ y Since its magnitude in graphene 

is taken to be around 2.7 eV in most cases, the adopted 
energy unit is about 6.3 eV. 

If chemisorbed atoms play the role of defects, even 
a minimal impurity model should contain, in a general 
case, a possibility for the electron transfer from the host 
to some energy level that belongs to the adsorbed atom. 
This task is fulfilled by the well-known Fano-Anderson 
impurity model, 17 which was adopted in quite a few stud- 
ies devoted to adatoms on graphene ^ 10 ' 18 On the as- 
sumption that impurities are deposited without any spa- 
tial correlation, we arrive at the following impurity part 
of the Hamiltonian: 

)], (2) 

where dimensionless 7 is the hybridization between the 
adatom and the host, Eq is the bare energy of the adatom 
level, djj a and d na are the creation and annihilation op- 
erators at this level, the variable r\ na takes the value of 
1 with the probability c or the value of with the prob- 
ability 1 — c, and c is the impurity concentration. 

By eliminating wavefunction amplitudes at adatoms 
from the stationary Shrodinger equation for the operators 
(TT|) and @ , the problem can be reduced to a more simple 
form with 



^ ^ Vna ^na^'n < 
na 



e-e 



(3) 



A conclusion on the presence of a resonance state for an 
isolated impurity can be made based on the single-site 
T-matrix, 



t(s) 
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where go(e) is the diagonal element of the host Green's 
function, 



(5) 



taken in the site representation. The required diagonal 
element can be approximated in a vicinity of the Dirac 
point (e = 0) as follows: 



g (e) « 2eln(|e|) - i7r|e|, 



lei « 1. 



(G) 



Substituting this expression into Eq. (j4]) and taking into 
account that logarithm is a slowly varying function, we 
have: 



r(e) 



l7e/| 



£-£ e/ + i7r|7 e/ | 2 |e| 



(7) 



where the effective hybridization and the effective energy 
of the impurity level are introduced: 



iTT _ _ £o 
D ' e/ ~ D ' 



I7e/| 

D = l-2| 7 | 2 (ln|e r | + l). 



(8) 



It immediately follows from Eq. (J7]) that e e f = e r , i.e. 
corresponds to the energy of a possible single-impurity 
resonance. 

This resonance can be considered as a well-defined, 
when its damping is significantly less than the energy 
interval between e r and the closest van Hove singularity 
of the host spectrum. In the case under consideration, 
the role of such singularity is played by the Dirac point. 
Thus, the mentioned condition for the appearance of a 
well-defined resonance is reduced to the inequality: 



8 = n\ lef \ 2 < 1, 



(9) 



where 8 is the small parameter that will be used further 
on. It is evident, that the inequality © is automatically 
satisfied at |-y | -c 1. Thus, a well-defined resonance state 
is always present in a system with Dirac-like dispersion, 
when the impurity level is weakly bound and the reso- 
nance energy is located close to the Dirac point. For a 
relatively large hybridization constant, \j\ > 1, the in- 
equality (J9j) takes the form: 



2|ln|e r | + 1| 



<1, 



(10) 



which can be satisfied only at a close proximity of the res- 
onance energy to the Dirac point, and coincides with the 
corresponding condition for Lifshitz impurity centers^* 
One can expect that at [7] > 1 the spectrum rearrange- 
ment, for the most part, also proceeds according to the 
anomalous scenario, which already has been examined 
in detail for Lifshitz impuritiesJ^^ Therefore, we are 
focused on the opposite case of weakly bound impuri- 
ties (| 7 | <C 1) below. In other words, the hybridization 
parameter has to be smaller than the transfer integral 
magnitude in the host. 



III. COHERENT SPECTRUM 
REARRANGEMENT 

At low impurity concentrations, c <C 1, the spectrum 
rearrangement can be analyzed by means of the mod- 
ified propagator method. 19 Within this approach, the 
self-energy cr(e), which approximates the averaged over 
impurity distributions Green's function of a disordered 
system 



G = 



1 



H 



1 



e-a(e)-H ' 



is determined self-consistently: 

c| 7 | 2 



a(e) = 



£0 



h\ 2 9o[e-*(e)Y 



(11) 



(12) 



We would like to remind that this method does not 
work satisfactorily in all spectral intervals of a disordered 
system. — Due to increased scatterings on impurity clus- 
ters, concentration broadening areas are formed nearby 
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singular points of the spectrum. Inside these areas the 
approximation ([T2"j) is not valid and states turn out to 
be localized. Thus, in order to locate such areas and to 
determine their widths, it is often enough to employ the 
standard Ioffc-Regel criterion^ 



|e(e)| = \e- Rccr(e)| > - lma(e), 



(13) 



where e(e) is the renormalized energy. According to the 
approximation (JTTJ), it specifies the renormalized disper- 
sion relation in a disordered system inside those spectral 
domains, which are occupied by extended states: 



e{e) = e(fe), 



(14) 



where e(k) is the host dispersion. It is well known that 
the Hamiltonian ([T]) yields the linear dispersion for free 
charge carriers, when the wave vector is counted from 
one of the two inequivalent Dirac points in the Brillouin 



e(k) ps ±VFk, 



where a is the honeycomb lattice constant. 

Let's assume that the spectrum rearrangement has oc- 
curred already. After the rearrangement, the spectrum of 
the disordered system should undergo a cardinal change. 
We are going to pay attention only to states that re- 
main extended despite the presence of disorder and to 
make several estimations that should help to sketch up 
the overall structure of the rearranged spectrum. Due to 
the singular character of the impurity perturbation, the 
distortion of the spectrum is most pronounced around 
the resonance energy e r . 

Taking into account expressions ([6][8]), and (fl"2j) . the 
imaginary part of the self-energy can be approximated 
as follows: 



Imcr(e) ps — 



«r|7e/|*|e(e)| 



( £ - £r )2 + (7r | 7e/ |2 £ - (e)) 



(16) 



when the condition (|13[) holds, and, consequently, the 
corresponding states can be treated as extended ones. 
While remaining inside the domain of extended states, it 
is possible, regardless of the strong impurity scattering 
near the resonance energy, to go further and to omit the 
second term in the denominator of (II 6[) (justification of 
this approximation will be provided below). Then, the 
inequality (|13|) in a vicinity of e r can be rewritten as: 



C7r|7e/| 4 |e(g)| 



which turns to 



\s - e r \ > v^lTe/l = A r . 



(17) 



(18) 



Thus, the width of the concentration broadening area 
of the impurity resonance A r is increasing proportional 



to y/c, i.e. behaves like the width of the mobility gap, 
which opens under the anomalous spectrum rearrange- 
ment in graphene with Lifshitz impurities; — However, it 
is worth noticing that, in comparison with the case of Lif- 
shitz defects, the width A r contains the small parameter 
of the problem |7 e /| 2 as a multiplier. Similarly, within 
the same approximations, we have inside domains of ex- 
tended states that 



i(e) 



C\leff 

e — E r 



(19) 



The renormalized dispersion relation e(fc), which is de- 
picted in Fig. [TJ can be obtained from Eqs. (fl4|) and (fi"9|) 
as a solution of the quadratic equation, 

e(k) ps i{ £r + s(k) ± .0 £r -e(fc)] 2 +4c| 7e/ | 2 }. (20) 

Obviously, positions of nodal points, i.e. those points at 
which dispersion cones are touching each other by their 
vertices, are determined in the rearranged spectrum by 
the condition e(k) = 0. Thus, it follows from the disper- 
sion law (|2"0")) that there should be a second (impurity) 
nodal point with the energy 



^nod 



2^ 



4c|7, 



ef\ 



(21) 



where we also assumed that e r > without loosing gen- 
erality. At low impurity concentrations, 



C\lef\ 

e r 



(22) 



moreover, A no d <C e r - On the other hand, one can speak 
about the appearance of the extra nodal point in the 
spectrum only in the case, when it lies outside the con- 
centration broadening area of the impurity resonance, i.e. 



(23) 



In essence, the condition (1231) determines the critical con- 
centration of the coherent spectrum rearrangement, 



7T£„ 



(24) 



It is not difficult to see that this magnitude agrees with 
quick-and-dirty estimates for the critical concentration, 
which are based on the mutual spatial overlap of individ- 
ual impurity states. Namely, the period of spatial oscilla- 
tions of the host Green's function in graphene determines 
the characteristic radius of the single-impurity state: 
Timp ~ ae^ 1 . With increasing the impurity concentra- 
tion, the average distance between impurities f ~ ac -1 / 2 
is gradually decreasing. Both lengths become equal at 
c ~ e 2 , which is in accord with the relation (|2~4"j) j 12 ' 14 



As it follows from Eqs. (fT5j) and p^|) . the absolute 
value of the renormalized energy e(e) is of the order of 
\fcpK at the boundaries of the concentration broadening 
area located at the resonance energy, and, according to 
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It follows from Eq. (|2Tj) and the definition ([22)1 that at 





// 






\ \ & r 





~~ ^max ft 


V\ ^max 




c » c* 



lTe/| 



(26) 



(a)c < c* 



(b)c > c. 



the width of the anomalous dispersion domain A noc i be- 
comes of the order of v^lTe/li and thus scales differently 
with the impurity concentration. Since c csr /c» ~ 5, the 
small parameter, which was defined by Eq. (J9)), this char- 
acteristic concentration always exceeds the critical con- 
centration of the spectrum rearrangement. At the onset 
of the spectrum rearrangement, the quantity A r has a 
square root dependence on c, while A noc ; has a linear 
one, which ensures its leading growth. By contrast, both 
widths have a square root dependence on the impurity 
concentration at c > c*, while their ratio remains nearly 
constant: 



FIG. 1: (Color online) Sketches of the dispersion of charge 
carriers in graphene with weakly bound impurities after 
the cross-type spectrum rearrangement drawn according to 
Eq. (|20[) : (a) c csr < c < c, the Fermi velocity at the im- 
purity nodal point is significantly reduced; (b) c, < c < 1, 
the Fermi velocity at the impurity nodal point is comparable 
to the one in the host system. The host dispersion is shown 
by the dashed line. Branches with the anomalous dispersion 
are given by the thick (green) solid line. The concentration 
broadening area is marked by the (yellow) shade. 




the relation (IM1) . far exceeds e r in the rearranged spec- 
trum. The corresponding maximum value of the wave 
vector k max (see Fig. [T]) that is reached at the bound- 
aries of the concentration broadening area appears to be 
of the order of the inverse mean distance between impu- 
rities, as it would be expected: 



(25) 



By substituting the magnitude ycpn in place of e(e) in 
the denominator of Eq. (fT6|) , one can easily verify that the 
second term is small in comparison with the first one by 
virtue of the inequality (|I8[) . This justifies the omission 
of the second term in the denominator of Eq. (fTrJf for the 
extended states. 

When energy is crossing e r , the corresponding e(e) s» 
c\"f e f | 2 /(e— e r ) changes its sign. In this way a spectral do- 
main with anomalous dispersion develops between e r and 
£nod (see Fig. 1(a) I. With respect to the valence band 
(as it appears in the Figure) , this domain can be consid- 



ered as an impurity band filled with extended (itinerant) 
states. Because of the particle-hole symmetry of the host 
spectrum, this impurity band fits exactly inside the gap 
that opens under the cross-type spectrum rearrangement 
of the conduction band. At that, the impurity band for- 
mation occurs simultaneously for both Dirac cones in the 
Brillouin zone of graphene. 



c> c* 



(27) 



In addition, the width of the spectral domain with 
ano malou s dispersion A no d exceeds \e r \ at c > c* (see 
Fig. [1(b) |. 



The Fermi velocity if at the extra nodal point £ no d 
can be expressed through the Fermi velocity of the host 
Vf by means of Eq. ([2U]) : 



Vp 



1 di(k) 



dk 



de{k) 



fe=0 



de{k) 



e(fc)=0 
1 



C 

4— 



-1/2 



(28) 



where the plus sign was chosen for the square root term 
in Eq. (|2"U1) . It is easy to see from the ratio (|2"51) that the 
Fermi velocity Vf is reduced in c/c* times at c < c 4 in 
comparison with the one of the host. As a result, it can 
be difficult to observe the branch with the anomalous 
dispersion in real experiments. However, this velocity 
is gradually approaching up/2 for c 3> c*, which elimi- 
nates indicated difficulties. Because the Fermi velocity 
near e noc i varies in several times, and the width of the 
anomalous dispersion domain varies with the impurity 
concentration according to a different law, the character- 
istic concentration c* separates two dissimilar regimes of 
the coherent spectrum rearrangement. 

The width A r is of the order of 6\e r \ at c ~ c csr , i.e. is 
substantially smaller than the interval between e r and the 
initial position of the host Dirac point. This relationship 
becomes reversed at 



c > c s 



Alef\ 4 



C-csr 



(29) 



Therefore, the asymmetry, which is present in the spec- 
trum due to the nonzero energy of the impurity resonance 
state, is effectively smeared out, when the impurity con- 
centration is that high. 
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In the first approximation, as it follows from Eqs. (fT2|) 
and (|22p . the impurities under consideration behave in 
a close vicinity of e no d at c csr < c < c* like Lifshitz 
defects with the on-site impurity perturbation vl ~ s r /c. 
According to the results obtained in Ref. [H, the critical 
concentration of the anomalous spectrum rearrangement 
for such defects can be roughly estimated as c asr ~ 1 jv\ . 
Thus, the condition c < c asr ~ c 2 /£ 2 should be satisfied 
for the system to remain in the unrearranged regime. 
This condition reduces to c > c csr , which is, for certain, 
fulfilled. It was demonstrated in Ref. [TU that in such a 
case the concentration broadening area of the Dirac point 
is exponentially small as compared to the bandwidth. 
Thus, its presence can be, to a known extent, overlooked. 

Similarly, the effective Lifshitz impurity perturbation 
in the vicinity of e no d is Vl ~ \lef\l \fc at c > c», and 
the exponential smallness of the corresponding concen- 
tration broadening area is guaranteed by the inequality 
c < c/|7 e /| 2 , which is fulfilled on account of the condition 
©. 

The host Dirac point, which is always present in 
the spectrum, is also shifted from its initial position 
due to the effect of impurities. As it directly follows 
from Eq. (l20l) . the magnitude of this shift is about 
— c|7 e /| 2 /|e r | at low impurity concentrations, which 
agrees with the width of the anomalous dispersion do- 
main at c > c csr . At that, this shift occurs in a direction 
opposite to the concentration displacement of e no d from 
e r . 

For c < c*, the effective magnitude of the impurity 
perturbation vl ~ — |7e/| 2 /l e i-l in a vicinity of the shifted 
host Dirac point, and the corresponding critical concen- 
tration c asr ~ £r/l7e/| 4 - The inequality c* < c aS r is 
again met in view of the condition ©. Thus, the width 
of the concentration broadening area remains once more 
exponentially small. At c > c*, the Dirac point shift is 
~ \/c|7 e /| by the absolute value, and the analysis of its 
concentration broadening can be performed in much the 
same way as it was done above for the extra nodal point 



vicinity. 



IV. CONCLUSION 

To summarize, with increasing the concentration of 
weakly bound impurity centers, which yield well-defined 
resonance states, the electron spectrum of graphene un- 
dergoes the coherent cross-type spectrum rearrangement. 
In contrast to the already known cases of the coherent 
spectrum rearrangement, it manifests the impurity nodal 
point in the spectrum and the impurity band with the 
anomalous dispersion. In addition, it features a concen- 
tration broadening area, or a mobility gap, around the 
impurity resonance energy, and thus, when the Fermi 
level position in the system is controlled by, for exam- 
ple, the gate voltage, it should be possible to observe a 
metal-insulator transition on Fermi level's entering the 
mobility gap, and then a re-entrant insulator-metal tran- 
sition on its leaving. For Lifshitz defects, the mobility 
gap is narrower than the bandwidth in ~ yfc times. In 
contrast, the mobility gap width for weakly bound de- 
fects A r contains additional small parameter as a multi- 
plier (see the inequality (O). Therefore, weakly bound 
impurities are more favorable for the experimental study 
of the re-entrant insulator-metal transition in graphene. 
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